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Abstract 



We extend Green's function method developed by Stewart and Gong to calculate the power spectrum 
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1 Introduction 



It is generally believed that the curvature perturbation produced during inflation jlj is the origin of the 
inhomogeneities necessary for large scale structure formation. The power spectrum of these perturbations 
is predicted to be approximately scale invariant by most of the inflation scenarios 2 , and this is confirmed 
by many recent cosmic microwave background observations and galaxy surveys |3j. The observation 
results from WMAP 0] , SDSS and other experiments give more accurate measurement of the power 
spectrum and the spectral index than ever before. Thus, it is very important to calculate the power 
spectrum precisely to match the observation results fully and to distinguish different inflation models. 

Many works have been done on the power spectrum for the density perturbations produced during 
inflation |^. For example, the power spectrum and spectral index of density fluctuation are calculated up 
to the first order of slow-roll parameters in (Jj ■ In jH| Stewart and Gong further set up a formalism that 
can be used to calculate the power spectrum of the curvature perturbations up to arbitrary order in the 
slow- roll expansion, and calculate the power spectrum and spectral index up to second-order corrections. 
Recently, Myung et al.Q have calculated the second-order corrections in the brane inflation model and 
in the noncommutative space-time inflation model, respectively. 

However, we note that in most calculations made before, the sound speed is assumed to be 
constant. Of cause, it is true in the inflation models derived by a scalar field of canonical kinetic term, 
where = 1. But, as we will see in the following, the sound speed is time-dependent in the general 
case. For instance, in the case of tachyon inflation (see for example |12lll8j and references therein), which 
has received a lot of attention in recent years, = 1 — |e, where e = —H/H^ is a first-order slow-roll 
parameter and H is the Hubble parameter, a dot denotes the derivative with respect to cosmic time t. 
If one only considers approximations up to the first-order corrections, because e is a second-order small 
parameter, e and then can be treated as constant approximately. But when one attempts to calculate 
second-order corrections, one has to consider the time-dependence of these quantities. As a result, the 
Mukhanov equation governing quantum fluctuations produced during inflation becomes more difficult to 
solve than the case with a constant sound speed. 

In this paper, we have extended Green's function method developed by Stewart and Gong [5] and 
obtained the second-order corrections in the slow-roll expansion to the power spectrum in the case with 
a time-dependent sound speed. In particular, the case of tachyon inflation is included. Note that the 
formalism can be used to calculate the power spectrum of the curvature perturbation up to arbitrary 
order in the slow-roll expansion. We use the units h = c = SttG = 1 throughout this paper. 

2 Preliminaries 
2.1 Basics 

Our starting point is the following effective action during inflation 




(2.1) 



where g is the determinant of the metric, R is the Ricci scalar, tp denotes a scalar field and 



(2.2) 



1 



is its kinetic term. The Lagrangian p plays the role of pressure and the energy density is given by 

p^2Xp,,-p. (2.3) 

Here denotes the partial derivative of p with respect to X. In this paper, we consider a flat universe 
for simplicity, namely K = 0. The background metric is 

ds^ = dt^ - a^{t)dx = a^ir]) {dr]^ - dx^) , (2.4) 

where rj is the conformal time defined by drj — dt/a. The action for scalar perturbation can be written 
down as ^Ol 

1 /■ r . „ . 7" .1 

drjd^x, (2.5) 



w + G„ wV w H w 



z 

where a prime denotes the derivative with respect to the conformal time ?/, the canonical quantization 
variable w = C is the curvature perturbation, and z is defined by 

(«) 

and the sound speed by 

= (2.7) 

P,x 

Using the Fourier transform of w, from (|2.5|l one can obtain the equation of motion for mode Wk^ i.e. the 
well-known Mukhanov equation jlO| 

< + (^C^e - Wk = 0, (2.8) 

where k is the wave- number for the mode Wk- Furthermore, from Eq. H2.6() and the dynamical equation 
of the Hubble parameter 

H^-\{p + p), (2.9) 
we can express z by the slow-roll parameter e = —H/H^, which is a first-order slow-roll parameter, as 

z = V2^ei/2. (2.10) 

2.2 Sound speed 

Obviously, the sound speed Cg defined by Eq. (|2.7|l is generally time-dependant for a general Lagrangian 
p{(f, X). In the slow-roll phase where X is a small quantity, compared to the potential of the scalar field 
ip. As a result, a generic Lagrangian is expected to admit an expansion of the form 

p{ip, X) = V{ip) + D{if)X + B{if)X^ + (2.11) 

From Eqs. H2.11|l . (|2.3() and (|2.7() . one has the sound speed 



ABX „, (BX\ 
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Ci^l-^+2A^—j , (2.12) 

up to the term of X squared. Substituting Eqs. H2.11|l and H2.3|l into Eq. Ij2.9|l . we get 

2BX^ + DX-H^e = 0. (2.13) 
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This can be considered as a quadratic equation of X and the solution satisfying the convergence condition 
DX > BX^ is 

up to the second-order correction of slow-roll parameter. Then, one can expand the sound speed in terms 
of the slow-roll parameter as 

Because the coefficients B and D are dependent of (/?, and in turn dependent of time, it is still very 
difficult to deal with Eq. I|2.8(l . From now on, we restrict ourselves to a class of special case where the 
sound speed has the form 

d = 1- cie- C2e^. (2.16) 

Here ci and C2 are two constants of order of unity. This form of sound speed includes some well-known 
examples. For instance, in the case of canonical scalar field one has ci = C2 = while ci = 2/3, C2 = 
for the case of tachyon inflation ^1 . 

In fact, by using the formalism developed in section 3, we can consider a class of more general sound 
speed as 

— I — Ciem — C2€peq, 

where is defined by Eq. H2.17|l and i = m, p, q are some positive integers. We will give a simple 
discussion on this in the appendix. But here, we will concentrate ourselves on the case with the sound 
speed given by Eq. H2.16(l . 

2.3 Slow-roll parameters 

In this paper, we introduce a set of slow-roll parameters expressed by Hubble parameter and its derivatives 
with respect to time as 

_ (-1)" /(P-^H 
namely, one has 

ei=e^--?r, 62 = -^. £3 = ... (2.18) 

evaluated at sound horizon crossing 10 , i.e. aH — Cgk. Note that, different from other slow-roll param- 
eters used for example in [HIE], these parameters defined in (|2.17|l are all first-order small parameters, 
namely satisfying |e„| < ^ for some small perturbation parameter ^. The higher order slow-roll parameters 
are products of the first-order parameters e„. In addition, it is easy to show 

en = e„ [ei + (-l)"+\e„ + e„+i)] , (2.19) 

which will be used extensively below. 

Note that for the case of canonical scalar field, the slow-roll parameters used in Ref. [5| are 



(2.17) 




and these have the following relations to the slow-roll parameters defined in (|2.17|l 

1 11 3 3 1 

ei = ei, ^1 = 2^2, 52 = --el- -^2^-3, h = ^4. + -ele^ - -e2e'iei, ... (2.20) 
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3 The calculations 



3.1 Slow-roll expansion and the formalism of calculations 

Our main task is to solve the Mukhanov equation Eq. (|2.8|) . Introducing new variables y 
X = — fc?7, we can recast the Mukhanov equation as 



0. 



2kwk and 



(3.1) 



Because is time-dependent, the problem is more complicated than the case with constant Cg. In order 
to obtain the solution of (|3.1|l up to second order corrections in the slow-roll expansion, we develop a 
formalism, in which the key point is to use Green's function method developed in Ref. [S] twice. The first 



step is to decompose the term 



where g{\a.x) is defined by 



z dx^ z 



into two parts and Eq. 1)3. l|l then can be rewritten as 



dx^ 



C 



y = A5(lna:)y, 



We can justify this decomposition as follows. From Eq. 1)2.10(1 we have 



:- = — + ( - I ( I - 2 ( - I ( ^ 



a 
a 



£1 



Cs 



Note that 
and 



— =a'H'{2-ei), 
a 



-krj : 



dt 



= -k 



da k 



(l-f ei -f 3ef + eiea) 



(3.2) 
(3.3) 

(3.4) 

(3.5) 
(3.6) 



Thus, for any C^, the leading term of 



lfz_ 
dx^ 



which comes from — . Therefore, the decomposition 



made above is always valid. Furthermore, we can expand 



g{lnx) = 



gn+i 



(Inx)" , 



(3.7) 



n=0 



where gn is of order n in the slow-roll expansion, namely |(/„| < This expansion is useful for 

exp(-l/^) < a; < exp(l/^). From Eqs. H3.3(l - (|3.7(l . H2.16|l and H2.19|l . up to second-order corrections, it 
is easy to find 



92 



dg 



din a; 



9 15 ^2 

-12el- yeie2 + 7^62 



262.3 



aH=C,k 



gi = g\x=\ 



3 3 11' 

6ei -t- (20 + 3ci)e2 + -62 + (9 -^c{)ex^2 - -^4 - ^"^^3 



(3.8) 
(3.9) 



aH=Csk 



Note that here In a; \ aH=Csk is a first-order small quantity in the slow-roll expansion. 

The boundary condition of Eq. 1)3.1(1 or Eq. ((3.21) required to get the power spectrum is 

lim yo{x) as a; — > 00, 
2kAkZ as a; ^ 0, 



(3.10) 
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where yo (x) is a solution of equation 



''"y ■ ^C!~\]y = 0. (3.11) 



Note that is time-dependent, and then x-dependant. To solve Eq. H3.11|l . it is necessary to express 
as a function of x. Since ei = —H/H^ is time-dependant and therefore is a function of x, too. From 
Eqs. (|2.19|) and H3.6|l . up to second-order, we have 

Here ei (2ei -|-e2) can be treated as a constant approximately because its derivative with respect to time is a 
third-order small quantity, which can be ignored. Therefore ei should be of the form of const. + const. In a;, 
up to second-order corrections. In fact, we can expand ei as 

oo 

^i = E^(l^^)"' (3-12) 

n=0 

where s„ is of order n in the slow-roll expansion, namely |s„| < This expansion is useful for 

exp(— 1/^) X exp(l/^). It is then easy to obtain 



32 



d\nx 



~ -ei(2ei-f e2)Lff=C.fc> (3-13) 

x=l 

•Si = - eilaff=c.fc- (3-14) 
Substituting Eqs. H3.12|l - H3.14|l into Eq. (|2.16|l . up to second-order, we can recast the sound speed as 

= Cf, -cisalnx, (3.15) 

in which 

C^^ = 1 - cisi - C2S^ (3.16) 
Then, Eq. H3.11|l can be approximately rewritten as 



^2.2 ' \ Csc-^jy = {ciS2\i^x)y. (3.17) 

Using Green's function approach. Equation H3.17|l can be perturbatively solved order by order. The 
solution of Eq. (|irTT|l or l|TT7|l is 

ici 

yo{x) = yo{x) + -TT^ du S2\nuyQ{u)[ya(u)yo{x) - yo{x)yQ{u)], (3.18) 



2Csc Jx 



where 



yo(^)= (l + 77^)e^^-^ (3.19) 



+ (cl-\]yo- 0. (3.20) 



satisfying 

dx'^ 

Ignoring the 4th- and higher order corrections, from H3.18|l we get 

yo{x) ^ yoix) + —- / du S2lnuyo{u)[yQ{u)yo{x) -yQ{x)yoiu)]=yo{x) +ysix). (3.21) 

^L-sc Jx 
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Now, we are able to solve Eq. (|3.2|1 by using Green's function approach once again. The solution is 

POO 



y{x) ~ ya{x) + 
- yo(a;) 



du\g{\nu)y{u)[yl{u)yn{x) - y^{x)yo{u)] 
du-^g{\nu)y{u)[yo{u)yo{x) - 2/o(x)yo(w)], 



(3.22) 



in which we have ignored the third and higher order corrections. Using Eq. 13.22|l twice, we have the 
perturbative solution up to second-order corrections 



y{x) = yoix) + ys{x) + yi{x) + y2i{x) + 2/22(2;), 
where ys{x) comes from Eq. (I3.21() . which is a second-order correction to y{x), and given by 
i ci 



ys{x) 



i C1S2 



dus2 \nuya{u)[yQ{u)yQ{x) ~ j/o(x)yo(u)] 



yo{x) 



2Cs. 



-1 — xhi x + X + —r^ (x ^ \nx — 1 + X ^) 



3i 
2Csc 



2/1 (x) is a first-order correction to y{x 

igi 



yi{x) 



2Cs 
1 



-351 



-yl{x 

r duu-^e^'^^^^ - (a;-ie2'^-"(lnx + l) -e 

X ^ sc 



duu '^yo{u)[yQ{u)yo{x) - yQ{x)yo{u)] 
^lLj:-^e'C.^-+y*oix) H duu-'e^'^-^ 

^sc Jx 



2Csc 



and the second-order corrections y2i(x) and 2/22(2;) are, respectively, given by 



2/21 (a;) = 



and 



2/22(2 



2Csc J 


1*00 

du 

X 




2 
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3Csc 


+iyoix) / 




J X 


ig2 - 

2Csc J 


du 

X 



duu '^yi{u)[y^{u)yo{x) - 2/o(x)2/o(u)] 



duu-^e'^'^''"" 



ig2 
3 



duu \a.uyo{u)[yl(u)ya{x) - yl{x)yQ{u)] 



^x-' (I + 2\nx\ e^^-- + ^y*{x) / duu-'e^^<^-- 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



(3.27) 



+iijl{x) / duu Mnue^^'^^^" 

J X 

Note that in writing 2/22(2;) we have ignored a term 

^ 922/0(2;) 3 d^lS] (lnx + - 

since it is a third-order correction. In addition, the perturbative solution Eq. 13.23(1 should be accurate 
for exp(l/0 > x > exp(-l/^). 



3.2 The power spectrum 

Having obtained the solution y{x) of the Mukhanov equation (|3.2|l . we are in a position to calculate the 
power spectrum of density fluctuation during inflation. Taking the limit x —* while keeping ^ln(f/a;) 
fixed and small, after some tedious calculations, we obtain the asymptotic forms for yo{x), ys{x), yi{x), 
2/21 (a;) and 2/22(2;). They are 

-1 



3jCiS2 _i 



2i 



4 

3 ~ 3a 
2/1 (a;) 



2iC,. 



2ie 



*5i -1 

X 



3 a 



3 Csc 

a + — -InC,, 



■a-2 + Sz(2iC,c)-y-lnC, 



»gi 

3 Csr. 



X """Inx, 



(3.28) 
(3.29) 
(3.30) 



y2iix) 



18 Csc 



a'-^a-4+^ + (lnC,e)2 + 



2a I In Csc + ITT a 



In a, 



9 a, 



1 1 ^ 

a---lnC., + y 



a; ^ Inx 



18 C„ 



x-^{\nxf 



(3.31) 



and 



V22{x) 



6 a, 



+ -a - — + iTT a + - - In Csc - 2 a + - In Csc + (InC^c)' 



*.92 -1, 

X mx 



9C 



V2 

ig2 
QCsc 



X~^{\Tixf, 



(3.32) 



respectively, where q; = 2 — ln2 — 7, 7~ 0.5772156649 is the Euler-Mascheroni constant, and Ei{z) 
—dx is the exponential integral function. Especially, Ei(2i) ~ 0.4229808288 + i 3.1762093036. 



The exact asymptotic form of y{x) in the limit a; ^ is given by Eq. (|3.1UI) . By comparing Ij^.tOfl 
with the sum of asymptotic forms given by (|3.28() - (|3.32fl . we can obtain the coefficient Ak- To this end, 
we need to expand z first, 

fn 



XZ 



5:^(ln.)", 



(3.33) 



where fn/fo is of order n in the slow-roll expansion. This expansion is useful for exp(— 1/^) ^ x ^ 
exp(l/^). By using Eqs. (|2.10() . H3.6|l . (|2.16() and H2.19|l . up to second-order, we obtain 



/2 



fi 



d{xz) 



din a; 
k I ei 

h' 



d^xz) 




(d Inx)^ 


x = \ 


k fei f 


hV' 


2\ 



k fei 

V 2" 



16ei + 96162 - 2*^2 - '^2e3 



aH=Csk 



77a/ ~ ( ~46i - 62 - 126^ - 4ci6i - 46162 - -C16162 



2 + (2 + ci)6i + ( 10 + 3ci + ^c^ + C2 ) 6^ + ( 3 + ici ) 6162 



aH=C,k 



y{x) V2kAkfox^^ + V2kAkfix-^ \nx + ]^V2kAkf2X-^{\nxf . 



(3.34) 

(3.35) 
(3.36) 



aH=C,k 

Then, up to second-order corrections, the asymptotic form H3.10|l for y{x) in the limit a; ^ can be 
expressed as 



(3.37) 
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Collecting Eqs. (|3.28|I - H,'132|I together also gives an asymptotic form for y{x) up to second-order correc- 
tions. Note that 



In a 



0(0+ o{e)- 



We can therefore throw away the third and higher order corrections in Eqs. (|3.28(l - (|3.32|l coming from 
Csc- Comparing the result with Eq. H3.37|l . the coefficient of will give Ak up to second-order correc- 
tions. The coefficients of x^^ lux and a::^^(lna;)^ simply give the consistent asymptotic behavior, namely 
proportional to z. We finally arrive at 



2k fo [a 
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1 _^ CiS2 



-2 - 4i + 3a + /3i + il32 



3i7r 



InC, 



+ 



18 



a - -a - 4 



TT 

T 



92 

' 6 



(3.38) 



where /3i + if32 = 2ie^' + 3Ei{2i) = -0.5496523673 z 8.6963342377. 
The power spectrum is defined by 









z 



27r2 



(3.39) 



Substituting Eqs. igSHll, (|XTBI) . (|^ . ifXT^ and (|XTH) into the above formula, we obtain 



Vsik) 



i/2 



1 + (4a - 2)ei + aea + 



7n^ 

4^2 - 23 — -h (2 - a - /3i)ci 



4 



^a^ -I- a - 11 



297r2 1 ^ ' 

^+2(2-a-A)c, 



£162 



^ + 12'^ 



(3.40) 



where the right hand side should be evaluated at aH — Cgk. The spectral index is defined by 



d\nk 



It is easy to get the result in present case. 



(3.41) 



ns{k) = 1 - 4ei - 62 + 8(a - l)e^ - ae2 + (5a - 3)eie2 - Q!e2e3 

7tt^ 

-4a^ + lOa - 27 + — + (3 - a - /?i)ci 



-31a^ + 60a- 172- 



1997r2 
12 



2 - -a^ - 



£2^3 



(20- 7a- 7/3i)ci 
^2 1 



e?e2 



TT 



7a^ 



?a + 22 - 



3l7r2 



12 



6 2 
{a + /3i- 2)ci 



Cl £163 



eie2e3 + 2 I 



12 



e2e3e4, (3.42) 



where also the right hand side should be evaluated at aH = Cgk. 
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4 Some special cases 



4.1 The case of canonical scalar field 

In this case, — 1, namely ci = C2 = 0. Substituting into Eqs. (|3.4U|) and H3.42|l . we will obtain the 
power spectrum and spectral index for the case of canonical scalar field. It is interesting to compare our 
results with the one obtained in Ref . . Having considered the relations between the slow- roll parameters 
in Ref. |H] and in this paper given Eq. H2.2()|l . and for the case of canonical scalar field, 

it is easy to see our results are completely the same as those of Ref. [S] . This gives a challenging check 
of our formalism. 

4.2 The case of tachyon infiation 

In this case C's = l-T^ = l--ei, namely ci = - and C2 = 0. Substituting these into Eqs. H3.40|l 

_____ o o 

and H3.42|l . we obtain the power spectrum and spectral index for tachyon case. Note that the first order 

corrections have been calculated, for example, in 13 . The result of the second-order corrections to the 

power spectrum and spectral index for the case of tachyon inflation is new. This is an concrete example 

that the sound speed is time-dependant. When one attempts to calculate the second and higher order 

corrections, can not be taken as a constant approximately. Thus, the formalism presented in this 

paper is applicable in this case. 

4.3 The first-order corrections 

When one considers only first-order corrections, the slow- roll parameters can be regarded approximately 
as some constants. From Eqs. (|3.40f) and 13.42|l . for any = 1 — ciei — C2ef , we have 

7',(fc)~-^[l + (4a-2)ei+ae2], (4.1) 

and 

ns{k) ~ 1 - 4ei - £2 + 8{a ~ l)ej - ael + {5a - 3)eie2 - ae2e3, (4.2) 

where again the right hand side should be evaluated at aH = Cgk. Note that in these expressions, 
the coefficients ci and C2 do not appear. This implies that up to the first order corrections, the time- 
dependence of sound speed does not have any effect on the power spectrum of density perturbation and 
spectral index. 

In summary, we have extended Green's function method developed by Stewart and Gong to calculate 
the power spectrum and its spectral index of density fluctuation produced during inflation in the case 
with a time-dependent sound speed, up to second order corrections in the slow-roll expansion. The result 
for the tachyon inflation is included as a special case. We have noted that up to the first order corrections, 
there are no effects of the time-dependence of sound speed on the power spectrum and spectral index. 
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5 Appendix 



As mentioned at the end of section 2.2, we can deal with a class of more general sound speed as 

= 1 - cie„ - C2epe,, (5.1) 

where m, p and q are some positive integers. Similar to the case of ei, from Eqs. H2.19|l and H8.6|l . we 
have 

^ ~ -x-h^ [ei + {-ir+He„, + e„+i)] (5.2) 

for any £„. And [ei + (— l)'"^^(em + e^+i)] can be treated as a constant approximately because its 
derivative with respect to time is a third-order small quantity which can be ignored. Thus integrating 
Eq. E2l yields 

Em = Smi + Sm2 InX, (5.3) 

up to second-order, where 

Sml — £m\aH=C^k, S„i2 — — Em [ci + (— l)™"*""^ (Cm + Cm+l )] | ^^^(-^^ ■ (5.4) 

We then can recast the sound speed as 

= Csc - ciSm2 Inx, (5.5) 

up to second-order, where 

C^c = 1 - CiSmi - C2SpiSqi. (5.6) 

The formalism developed in section 3 is still valid for this because the calculations keep the forms of 

Csc, Sm 5n and /„ unchanged. The only difference is to replace the S2 in Ak by Sm2- Because the variable z 

z" 

is dependent of Cs as in Eq. 12.1()ll . we must calculate the slow-roll expansions of z and ^(Ina;) = —tt 2 

accordingly. Once having the new /o, gi and 52, and using the new C^^ given in Eq. H5.6|l and Smi, Sm2 in 
Eq. (|5.4() . we then can get the corresponding Ak- Thus the final power spectrum and spectral index are 
in hand. The final results are very involved and their expressions are long in length, we do not present 
them here. 
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